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The luminosity distance-redshift relation is one of the fundamental tools of modern cosmology. We 
compute the luminosity distance-redshift relation in a perturbed fiat matter-dominated Universe, 
taking into account the presence of cosmological inhomogeneities up to second order in perturbation 
theory. Cosmological observations implementing the luminosity distance-redshift relation tell us 
that the Universe is presently undergoing a phase of accelerated expansion. This seems to call for a 
mysterious Dark Energy component with negative pressure. Our findings suggest that the need of 
a Dark Energy fluid may be challenged once a realistic inhomogeneous Universe is considered and 
that an accelerated expansion may be consistent with a matter-dominated Universe. 

PACS numbers: 98.80.Cq, DFPD 04/A/27 



I. INTRODUCTION 

One of the fundamental relations in cosmology is the one 
expressing the luminosity distance dr, of a cosmological 
source in terms of its redshift z. In recent years this re- 
lation has been exploited to measure the present value 
of the expansion rate, Hubble's constant, with increasing 
accuracy. With the exploration of the Universe at red- 
shifts of order unity, we now have information about the 
time evolution of the expansion rate . A most surpris- 
ing result is that the time evolution of the expansion rate 
does not seem to be described by a matter-dominated 
Friedmann-Robertson- Walker (FRW) homogeneous cos- 
mological model of the Universe. The usual explanation 
for the discrepancy is that there is a new component of 
the energy density of the Universe, known as Dark En- 
ergy, that determines the recent evolution of the expan- 
sion rate. Of course all indications for Dark Energy are 
indirect; from cosmological data we only infer that the 
Universe is presently undergoing a phase of acceleration. 

Since the luminosity distance-redshift relation is of 
such fundamental importance, we must understand any 
possible effects that would result in a relation different 
from the FRW prediction for a homogeneous Universe. 
The computation of the luminosity (or angular diame- 
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ter) distance in a locally non homogeneous universe was 
first addressed by Zeldovich 0). Dyer and Roeder Q 
adopted the so-called empty beam approximation to de- 
rive an equation for the angular diameter distance. Our 
techinique instead relies on the traditional perturbative 
approach, recently used in Refs. |f|. 

In this paper we study the change in the luminos- 
ity distance-redshift relation due to cosmological per- 
turbations present in a background matter-dominated 
Universe with zero cosmological constant. In particu- 
lar, we provide the generic expression for the luminosity 
distance-redshift relation at second order in perturba- 
tion theory, for any value of the physical redshift and 
any direction of observation. We only consider modi- 
fications to the luminosity distance-redshift relation of 
a matter-dominated Universe, although our results can 
be easily extended to a Universe containing a mixture 
of matter and other fluids. Our findings may be used, 
for instance, to estimate the influence of lensing on the 
brightness of supernovae sources. What our second-order 
procedure does not account for is another class of terms, 
where a first-order small deflection in the light ray leads 
the geodesic to deflect to a region where the perturba- 
tions are sizeable relative to the neighbourhood of the 
geodesic in the unperturbed space-time. 

It is well-known that the luminosity distance-redshift 
relation allows to extract the theoretical predictions for 
the local Hubble rate !H and deceleration parameter q 
upon expanding around the observer point at z — 0. In a 
perturbed Universe, however, both the Hubble constant 
and the deceleration parameter lose their deterministic 
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nature; one has to consider the statistical nature of the 
vacuum fluctuations from which the present-day gravita- 
tional potential is originated. Therefore, the theoretical 
predictions for the expected values of the cosmological 
parameters are accompanied by a nonvanishing cosmic 
variance implying an intrinsic theoretical error. When 
calculating the variance of the deceleration parameter, 
we uncover an interesting infrared effect. We observe, at 
second order in perturbation theory, a large contribution 
to the variance from the cosmological perturbations with 
the largest wavelengths. If inflation is the origin of the 
cosmological perturbations pj , the beginning of inflation 
provides a cut-off to the infrared modes of the fluctua- 
tions. Our results suggest that if the super-Hubble modes 
have physical significance, they could yield a most im- 
portant modification to the deceleration parameter. One 
might even speculate that a complete treatment of the 
effect could obviate the need for the Dark Energy as- 
sumption. Indeed, despite the fact that the homogeneous 
FRW model for the matter-dominated Universe and no 
cosmological constant predicts that the latter may not 
decelerate, i.e. qo = | > 0, because of the large vari- 
ance, the true locally-determined value of the decelera- 
tion parameter has a non-zero probability of being less 
than zero. Put it differently, the theoretical prediction in 
a perturbed matter-dominated model is not a single well- 
defined curve in the luminosity-distance redshift plane, 
but instead is represented by a finite confidence region 
whose size is determined by the cosmic variance. 

The paper is organized as follows. In the next section 
we obtain the generic luminosity distance-redshift rela- 
tion valid at any order in perturbation theory around the 
homogeneous FRW background. In Section III we ex- 
press the luminosity distance-redshift relation in terms 
of the metric perturbations expanded up to second or- 
der. Section IV is devoted to the evaluation of the mean 
and the variance of the cosmological observables and to 
the discussion of their implications. Finally, Section IV 
present our conclusions and comments. All the technical 
details are included in Appendices A and B. 



II. THE GENERIC LUMINOSITY 
DISTANCE-REDSHIFT RELATION 

The goal of this section is to provide the reader with the 
generic formulation of the luminosity distance-redshift 
relation in a general setting. Our treatment closely fol- 
lows the one given in Ref. Q. In the next section, we 
will make use of this generic formulation to obtain the lu- 
minosity distance-redshift relation in a perturbed FRW 
Universe being the background matter-dominated Uni- 
verse with zero cosmological constant. 

In order to deal with the propagation of light from a 
given source to the observer, we make use of the confor- 
mal (Weyl) invariance of the electromagnetic field and 
define 



where x M = (77, x l ) are the space-time coordinates and 
the scale factor a is normalized to unity at the present 
conformal time 770 (a (770) = ao = !)• 

We will use a "~" to mark quantities calculated in a 
space-time with the metric g^ u ; quantities without """ 
will be calculated with the metric g^ v instead. Further- 
more, quantities with a "~" stand for quantities in a per- 
turbed Universe. 

Using the geometric optics approximation, the energy- 
momentum tensor of a photon emitted by a given source 
is 



T" v = A 2 k»k v 



(2) 



where A is the scalar amplitude of the wave and fc M is 
the photon four-momentum: 



dr. 1 * - - 
dv 
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(3) 



Here S is the phase (eikonal) of the wave and v is an 
affine parameter along the ray (i.e. along the photon 
trajectory) . The evolution of the parameter A is provided 
by the continuity equation V ' V T^ V = 0: we have 



(4) 



All these equations can be replaced by the corresponding 
ones in the metric g^ v . If we define 



d\ = a 2 dv 
we easily get 



dX 



(5) 
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dX 



dX 2 

1 
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af3 dX dX 
Aa6 , 9 = V M fc M 



(6) 



To obtain the equation describing the transport of the 
optical scalar 9 (expansion of the null congruence fc M ) 
along the ray, it is useful to recall the concept of null 
(Newman- Penrose) tetrad: it is made of two real (fc M and 
m M ) and two complex conjugate (t^ and i M ) 4-vectors, 
satisfying 



Pm^ = i% = 1 



.„ - m"m M = t% = A% = m% = . (7) 

Defining the projector 

3>£ = tH v + tuP =6£- k u m» - m v k» , (8) 

one can easily obtain the decomposition (the semicolon 
indicates covariant differentiation in the conformal space- 
time with metric g^) 



ds 2 = g^dx^dx" = a 2 (77) g^dx^dx" 



(1) 



(9) 
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with 



A^ v — k^pf ^ y^, — 25ft y ~^tfitu + ot^t v 



, (10) 

where 9 is real and a is complex. With such a decompo- 
sition we obtain 

k^k^ - ?1 ■ (ii) 



kl 2 = 



l 



2M 



(13) 



The quantity a is called shear of the null congruence. 
/From the Ricci identity 

k KV s -(y<-+8) = Ra^usk" 7 , (12) 

contracting the indices /i and u, we obtain (see for in- 
stance y|) 

dX - k 2 

In a similar fashion, the transport equation for the shear 
<7 can be derived by contracting (|12f> with t^t^k 5 and 
assuming that i M is parallel-transported along the ray; 
we get 

^ = -aO + C a0l/tl k a kH u tP, k v V v ¥ = , (14) 
dX 

where C a p vll is the Weyl tensor. 

To calculate the luminosity distance, we need the en- 
ergy flux per unit surface, I, measured by an observer 
with 4-velocity u M 



£=JKu (u a f°») (u s f s ") 



where 



(15) 



(16) 



If the source has physical radius R (which we will even- 
tually set to zero) and if we choose the affine parameter 
A on the ray connecting the observer and the source such 
that A = 0, A = A s , A = A s + AA S correspond to the 
observer today, to the surface of the source and to its 
center, the power & emitted by the source and measured 
by a comoving observer, is 

£ = AnR 2 i{X s ) 



The luminosity distance thus reads 



'£(X S 



R- 



i{0) ~A(0) 
where the redshift z is defined by 

w(A.) 



(1 + , 



l + z(X s ) = 



w(0) 



(17) 



(18) 



(19) 



To summarize, the evaluation of the luminosity distance 
goes through the following steps: one has to solve for 
the photon trajectory from which one may deduce the 
frequency u; then one is able to solve for the expansion 
parameter 6 and the shear a; finally one solves for the 
amplitude A. 



III. THE LUMINOSITY DISTANCE— REDSHIFT 
RELATION IN A PERTURBED FRW UNIVERSE 

The formalism we have summarized in the previous 
Section can now be applied to obtain the luminosity 
distance-redshift relation in a perturbed FRW Universe. 
In particular, we are interested in a perturbative treat- 
ment up to second order. This Section contains only the 
main steps of the calculation, a plethora of details can be 
found in Appendix A. 

First of all, we need to expand the wave four- vector of 
the photon reaching the observer at second order, 



* — (0) + ""(1) ^ K (2) 



(20) 



Similarly, we expand the photon trajectory as 

x*(\)=xfa(\)+x* 1) (\) + xfc(\). (21) 

Both quantities k^ and the (r — 0,1,2) are fully 
determined by the geodesic equation 



(kt, +ife£, +k? 



dX V (°) ' ^ ( 2 ) 



1 (l)a^ ft (0) /i (0) 



( r M h a 

\ L (2)a/3 ft (C 



(0)^(0) 



+ 2r (l)a/3 fc (0)^fl) + ^(l)a0,y k (O) k (O) X Jl)) ' ( 22 ) 

and by the initial conditions (given in Appendix A). We 
then need to perturb the expansion 6 of our null congru- 
ence 



7(0) 



1(2) 



(23) 



these functions are determined by Eq. i|13|l and by the 
boundary conditions. Furthermore, since the right-hand 
side of Eq. (|13fl contains the square of the shear, we 
need also to solve Eq. l(Ti|) up to first order, again with 
suitable boundary conditions. After the quantities 8, a 
and A are known, one is ready to calculate the luminosity 
distance-redshift relation. Sending the physical size R 
of the source (and therefore AA S ) to zero, we get the 
following expression 



d L = (l + z(A s ))a A 
1 




(1) (A, A s ) dX 



e {2) {x, x s ) dx - k° {1) (x s ) 



k° {2) (\s) + ^k° {1) (X s ) d (1) (X,X s )dX 



A, 



.(24) 



This is the generic expression for the luminosity distance- 
redshift relation at second order in perturbation theory 
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and for any value of the physical redshift z(X s ) and any 
given direction of observation. 

Since the physical observable is the redshift and not 
the affine parameter A s , we have to trade A s in terms of 
z in Eq. 124|) . In order to achieve this we may make 
use of the definition of the redshift in the unperturbed 
Universe 



where both X(i)ij and X(2)ij are traceless. 

Assuming that the source and the observer are both 
comoving with the fluid, one has (at any order in pertur- 
bation theory) w M = 5q . By solving Einstein's equations, 
the metric perturbations %j)i r > X(r)> r — 1, 2, can be ex- 
pressed in terms of the peculiar gravitational potential tp 



l + z 



«0 



a(r)o - A s ) ' 770 



(25) 



from which we deduce the relation between the physical 
redshift z and z 



1 + Z{K) 



ao {u^) (A s ) 
a( V (X 3 )) (u M fr*)(0) 

= ao(l-(k (1) (\ s ) + k? 2) (A s ))) 



a (V(o) (A s )) + a' (t? (0 ) (A s )) (t? (1) (A s ) + rj (2) (A s )) 



-a" (r? (0) (A s )) 77(1) (A s ) 2 



= (l + z)(l + T 1 (A s )+T 2 (A s )) , 



(26) 



where ij^ and ?7( 2 ) are respectively the first- and second- 
order expansion of the photon conformal time (see Ap- 
pendix A) and 

^(A)=(-^ (1 )-4))(A). (27) 



T 2 (A) 



o'fa(o)) , «"(»7(o)) 2 

^(2) " fc (2) ~ W^TT^l) 



+ 



We may therefore compute first the luminosity distance- 
redshift relation in terms of the parameter z and then 
as a function of the physical redshift z using Eqs. (j26|) . 
(|27|l and l|28|) to express z as a function of z. Yet, in the 
Appendix we describe another procedure, which allows 
to obtain an explicit expression of as a function of the 
physical redshift z. Both give rise to the same result. 

To compute all these quantities we have to expand the 
metric to second order. We adopt the comoving and syn- 
chronous gauge 1 



ds 2 



-drj 2 + "fijdx l dx^ 



lij = (1 - 2^(1) - V>(2))<% + X(i)ij + 2*(2)y > ( 29 ) 



1 The luminosity distance-redshift relation is a gauge-invariant 
concept as explicitly shown in Ref. Q and therefore the choice 
of the gauge is arbitrary. 



X(i)ij 



-tp- 



18 y ' 



SijV 2 tp 



at first order and 

, 50 2 5i] 2 k 77 

V>(2) = — g-V —54 ^' <P,k+ 



252 



10 



<P' ik <P,ik + 



(30) 
(31) 



(32) 



(2) /in 

19 
~3 



12,=.,, V 2 ,; + 4(VV)% 



19 fci r \ 10r/ 2 / 1 fe . \ 



(33) 



at second order 2 . We have disregarded both linear vec- 
tor modes (since they are not generated during inflation) 
and linear tensor modes (because their dynamical role is 
negligible). At second order tensor modes described by 
TTij do not enter in the computation of the Hubble rate 
and the deceleration parameter. Furthermore, the values 
of the second-order potentials have been computed with 
a proper match to the initial conditions set by single-field 
models of inflation [lG|. 3 

In Appendix A we will give explicit expressions for 
fc^(A) and 0( r ), r = 1,2 in terms of tp: in particular 
we will see that the luminosity distance d^, given by 
Eq. (|A69(1 . involves, besides the physical redshift z and 
the direction of observation, only the peculiar gravita- 
tional potential tp and its gradients along the background 
geodesic connecting the source to the observer and the 
present conformal time 770 = 2/Jfo. 



IV. RESULTS AND IMPLICATIONS 

We are now ready to extract from the findings of the 
previous Section the relevant cosmological parameters in 
a realistic perturbed Universe. To this purpose, let us 



2 The well-known residual gauge ambiguity of the synchronous and 
comoving gauge has been fixed as in Ref. 

3 Notice, however, that the initial conditions do not play a signifi- 
cant role, since physical observables like the Hubble constant and 
the deceleration parameter depend only upon the rate of change 
of cosmological inhomogeneities. 
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recall that in a generic unperturbed FRW Universe the 
following relation holds 



dL = W 



z+y(l-<7o) + Q(z 3 ) 



(34) 



where 3{q is the unperturbed Hubble constant and go 
is the unperturbed deceleration parameter. Therefore, 
an observer measuring a relation of the type dr, — 
A + H>z + Cz 2 + ■ ■ ■ would conclude that 23 and C provide 
a measure of the Hubble constant and deceleration pa- 
rameter at the present epoch. To make contact with this 



procedure, we may expand the generic formula l|24() for 
the luminosity distance (valid at any value of the redshift 
2) around z = 0. In this way, we may determine the value 
of the Hubble constant and the deceleration parameter 
that would be measured by an observer having at her/his 
disposal sources with redshifts z < 1 in a perturbed Uni- 
verse. Performing properly the angular average (• • • )n 
and comparing this expansion of Eq. (|24[1 with Eq. H34|) , 
we infer the expression for 3~Cq and for gb (see Appendix 
B for details): 



Jf 



2/1 i l / ■• y 

1 ~ K~ + 4^ 2 ) + + 12 ( X « J X(1)l 



IvV - — (V^) 2 + — ^V 2 ^ ( — 
18 ^ 108 1 27^ V V^o 



2 



+ ^ ( 2 V>(i) + V>( 2 ) + 4^(1)^(1) + I X(i)i?J 

^) f^w + i^( 2) + j (^[i)) 2 + ^(1)^1) + ^ (x ( i )y )' (*&)' + y 2 (xl)"xm 

^) (^ 1 )^ 1 )+i(x(i) ii )'(x;i ) )" 



1 + [ JLvV + — yvV- — (v^) 2 

1 18 * IT 108 v ^ 



2 



1 

30 



72 \ 2 23 
270 



2 



(35) 



(36) 



where "0(r)j X(r) ( r = 1, 2) and ip are evaluated at the 
observer's position xrm = and at the present time rjo. 
Eqs. H35|) an< i CUt are among the main results of our 
paper. The key point is that the values l|35|) and i|36|) 
of the Hubble constant and the deceleration parameter 
in a perturbed Universe are not deterministic. Indeed, 
we must consider the statistical nature of the vacuum 
fluctuations from which the present-day linear gravita- 
tional potential if is originated. This implies that the 
gravitational potential does not have well-defined values, 
but one can only define the probability of finding a given 
value at a given point in space. In a given realization 
of the perturbed Universe, the values of the cosmological 
parameters may change in regions of the Universe which 
are causally disconnected today. Therefore, it is unavoid- 
able that the theoretical predictions for the values of the 
cosmological parameters come with a nonvanishing cos- 
mic variance which implies an intrinsic theoretical error. 



4 Note that Eq. agrees with Eq. (39) of HI . 



It will turn out that the variance can be large and get 
its largest contribution from the infrared super-Hubble 
modes of the metric perturbations. 

Let us evaluate the expected uncertainty in the de- 
termination of the deceleration parameter, we treat the 
metric fluctuation ip as a Gaussian random variable with 
zero statistical mean over a volume which, in the infla- 
tionary picture, has a size much larger than the present- 
day Hubble radius. The variable ip takes random values 
over different "realizations" . We may express the vari- 
ance in terms of the matter power spectrum. The proce- 
dure would be to fix a spherical domain of large volume V 
surrounding the observer and containing the most distant 
sources of interest 5 . However, since the main contribu- 
tion to the variance comes from the longest wavelengths, 



5 We are interested in comparing the theoretical predictions of a 
perturbed Universe to the ones obtained in a unperturbed FRW 
homogeneous classical background. Since the quantities relative 
to the unperturbed model are recovered in the zero-momentum 
limit, the volume average has to be performed over a volume of 
size at least as large as the present-day Hubble radius. 
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this volume averaging is irrelevant 6 . We will express ip 
and its derivatives in terms of a Fourier integral, so 



d 3 k i 



d 3 k 



iki e l 



d 3 k 



k 2 n e lk ' s , 



etc. 



(37) 



(2tt) 3 

The Fourier components ipj: satisfy 
W k = , n^T 2 = {2^) 3 5 {3 \k 1 + k 2 ) P v (ki) (38) 

n 1 fk 3 n 3 ( Piu = (^) 6 {s (3) (ki + k 2 )s^(k 3 + k 4 ) 

x P v (ki)P v (k 3 ) + (cyclic terms)} , 



where (• • • ) denotes the statistical average and P v (k) is 
the power spectrum of the gravitational potential. We 
can express Pm(k) in terms of the matter power spectrum 
as 



p#)4« A2(u) 



k 7 



(39) 



where A 2 (A;, an) is the (dimensionless) power spectrum 
of the matter density fluctuations linearly extrapolated 
to the present time. We express the power spectrum 
A 2 (fc,an) in terms of the transfer function T 2 (k). For a 
Harrison-Zel'dovich spectrum, the power spectrum is 



T 2 (k), 



(40) 



where A is the dimensionless amplitude, ^4 = 1.9xl0 -5 . 

An analysis similar to the one performed in Ref. [Tl) 
shows that the biggest contribution to the variance of the 
deceleration parameter comes from the terms of the type 
ipV 2 tp whose variance is 

Var^VV] * QaX) 2 /^A 2 (fcx,a ) 
-A 2 (fc 2 ,a ) 



dk 2 A 2, 

k 5 



Var [ip] Var [V 2 ^] 



(41) 



Taking into account that T{k) — > 1 when k — > 0, we 
conclude that the variance of the deceleration parameter 
is sensitive also to the infrared modes of the Harrison- 
Zel'dovich power spectrum and is 



V Var [(g) n] 
<7o 



A 2 In - 



^■I IN 



10- 10 In 



«max 



>-MIN 



(42) 



As a technical remark, we stress that the low-pass filtering pro- 
cedure over the volume V should have been performed before 
Taylor expanding the luminosity distance— redshift relation 1241 
around the observer point at z = 0. This amounts to cutting off 
the spurious contribution of ultraviolet modes to the means and 
variances of physical observables which would appear performing 
the average after expanding in powers of the physical redshift. 



where /cmin is the infrared cut-off and &max is the ul- 
traviolet cut-off set by the averaging volume. We can 
take it to coincide with the horizon volume, that is 

^MAX = ■ 

Rather than a Harrison-Zel'dovich spectrum, if we as- 
sume a slightly red spectrum so that A 2 (k) cx fc 3 +™= with 
< (1 — n s ) <C 1, then the logarithmic term in Eq. I|42|) 
is replaced by (1 — n s ) _1 (/cMAx/fcMiN)^ 1 ~ ns - ) . Now this 
will give unit variance if 



In k M Ax/k 



MIN 



(45 + m(l - n s ))/(l - O . (43) 



The variance of the deceleration parameter is infrared 
sensitive. There are various ways to avoid such a sensi- 
tivity: 

• If the power spectrum is blue, i.e. n s > 1, the 
variance does not pick up a large contribution from 
the long wavelengths; 

• One may imagine to fix the physical infrared cut-off 
at a momentum scale corresponding to the present- 
day Hubble radius, but in the inflationary picture 
there is no strong motivation to do so since infla- 
tion is likely to have lasted for a period much longer 
than the minimum required number of e- folds ~ 60. 
After inflation, the size of the Universe which is rel- 
atively homogeneous is of the order of ~ Hj e N , 
where N is the total number of e-folds and Hi is the 
Hubble rate during inflation. As the length scale 

is much larger than the present-day Hubble size, 
it seems reasonable to sum up all the super-Hubble 
modes up to a length scale corresponding to £* at 
the end of inflation 7 . 

• The presence of an infrared sensitivity might be 
only an artifact of the perturbative expansion and 
it might disappear when dealing with a complete 
nonperturbative approach. To check this possibil- 
ity, a computation beyond second order is needed, 
even though we do not see any a priori reason why 
such cancellation should manifest itself at higher 
orders. 

Since our findings reveal the crucial role played by the 
infrared modes, a deeper understanding of their physi- 
cal significance is certainly desirable 8 . A large variance 



7 One should also keep in mind that on scales larger than t* the 
Universe may become extremely inhomogeneous due to quantum 
fluctuations produced during inflation and a considerable part of 
the physical volume of the entire Universe may remain forever in 
the inflationary phase 1 1 (I . 

8 For references dealing with the physical significance of the super- 
Hubble modes neglecting the gr adie nts see Refs. Il 71 and ll'-l . 
Recently, it is has been claimed lid that "the relativistic zero- 
pressure fluid perturbed to second order in a flat Friedmann back- 
ground coincides with the Newtonian result" and that "there are 
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of the deceleration parameter V Va ^ g ^" I is caused by 
the fact that the cosmological perturbations on super- 
Hubble scales are time-dependent when gradients are 
consistently taken into account. Indeed, it is easy to 
convince oneself that, were the infrared part of the per- 
turbations i/jm, X(r) ( r — 1)2) constant in time, they 
might be eliminated from the metric 129|) by a simple 
rescaling of the spatial coordinates which allows to re- 
main in the synchronous gauge. Such a freedom is lost if 
the inhomogeneities have a non-trivial time-dependence 
on super-Hubble scales. Consider, for instance, the first- 
order perturbation t/>m = + fg V 2 ip. It contributes 
to the metric <|29[) by a piece (l — ^§-<f) &ij which can be 
rescaled to &y by a transformation of the spatial coordi- 
nates. However, this change is not for free, it gives rise to 
the crucial piece proportional to ipV 2 tp at second order 
as it is easily realized by inspectingthe transformations 
of the perturbations listed in Ref. Q ■ The fact that the 
variance gets its largest contribution from a piece propor- 
tional to ip\/ 2 ip does not then come as a surprise, it man- 
ifests the impossibility of rescaling out the super-Hubble 
modes. An alternative way of dealing with super-Hubble 
modes can be found in Appendic C. 

In the following, we set the infrared cut-off to the value 
fixed by the beginning of inflation which, in turn, depends 
upon the total number of e-folds N of the inflationary 
period. 

The ratio (/cmax/^min) is predicted by the inflationary 
theory to be (we are adopting natural units) 



&MAX 



iMIN 



10 



-30 



Trr\ n 



(44) 



where Trh is the reheating temperature at the beginning 
of the radiation era after the end of inflation. 

For the variance to be of the order of the background 
value go = \, the perturbation spectrum would have to 
extend to a factor of exp(6 x 10 18 ) (10 188 e-folds!) times 
the present Hubble radius (for Trh — H\) in the case of 
the Harrison-Zel'dovich spectrum. However, if, for in- 
stance, n s = 0.94 on super-Hubble-radius scales, then a 
variance of order unity is obtained if the perturbation 
spectrum extends N ~ 700 e-folds beyond the Hubble 
radius. Since the present Hubble radius corresponds to 
a scale that crossed the Hubble radius about 60 e-folds 
before the end of inflation, if inflation lasted more than 
~ 700 e-folds with a super-Hubble-radius spectral index 
of n s = 0.94, then the effect of super-Hubble-radius per- 



no relativistic correction terms even near and beyond the horizon 
to the second-order perturbation" . However, this claim relies on 
disregarding second-order tensor mo des and is not justified, as 
implied by the discussion of Ref. f Il9l 1. where it was stressed the 
importance of second-order tensor modes for the correct recovery 
of the Newtonian limit from the full relativistic theory. 



turbations on the locally observed value of the decelera- 
tion parameter would be sizeable. 

What are the practical implications of our findings? 
What observations tell us is that the Universe is presently 
undergoing a phase of accelerated expansion, i.e. that 
the deceleration parameter is negative. Indeed, the un- 
expected faintness of high-redshift Type la supernovae 
(SNe la), as measured by two independent teams [l|, has 
been interpreted as evidence that the expansion of the 
Universe is accelerating. In an unperturbed FRW Uni- 
verse, the deceleration parameter is uniquely determined 
by the relative densities of the various fluids with their 
own equation of state 



1„ 3 r 



(45) 



where Slo is the present-day total energy density param- 
eter and fij are the relative contributions of the various 
components with equation of state Wi — Pi/ Pi (P and p 
are the pressure and energy density, respectively) . There- 
fore, the observation of a negative value of the decelera- 
tion parameter seems to call for the presence of a "Dark 
Energy" component with negative equation of state [l3| 
and abundance J1a ~ 0.7. 

The need of a mysterious Dark Energy fluid seems to be 
challenged once a realistic perturbed Universe is consid- 
ered. Our results show that the theoretical predictions 
for the local Hubble rate and the deceleration param- 
eter are affected by a cosmic variance whose size may 
be large, depending upon the value of the spectral in- 
dex and the overall duration of inflation. Suppose, as we 
have done so far, that the Universe is globally flat and 
matter-dominated with Qm ~ 0.3 and the background 
value go = h > 0. Because of the large variance, how- 
ever, the true locally-determined value of the deceler- 
ation parameter has non-zero probability of being less 
than zero! In other words, in a perturbed Universe, ac- 
celeration might not imply the existence of Dark Energy. 
Fig. 1 gives a qualitative sketch of the effect of the cosmic 
variance predicted in our paper as far as the magnitude- 
redshift relation is concerned: the theoretical prediction 
in a matter-dominated model is not a single well-defined 
curve, but instead is represented by a finite region whose 
size is determined by the cosmic variance. 



V. CONCLUSIONS AND COMMENTS 

Let us close with some comments. First of all, let us 
notice that in a perturbed Universe the theory is unable 
to provide the expected value of the Hubble constant 1)3 5|l 
and deceleration parameter (|36f) since they depend upon 
the "bare" unobservable Hubble constant JCq. There- 
fore, in order to compute the probability that a typical 
observer measures values of the Hubble constant and de- 
celeration parameter in agreement with the observations, 
a complete procedure would require marginalizing over 
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z 



Figure 1: The residual Hubble diagram with respect to an 
empty Universe for a collection of type la Supernovae taken 
from Ref. [20|. The apparent magnitude m and the absolute 
magnitude M are related to the luminosity distance by the 
relation (m — M) = — 5 + 51og 10 dk(pc). The red long-dashed 
line corresponds to go = 1/2 for a fiat matter-dominated Uni- 
verse (up to 0(z 3 ) terms); the blue dot-dashed lines corre- 
spond to the values q = qo ± 1, i.e. to the possible size of 
the cosmic variance band expected in a perturbed Universe; 
the orange short-dashed line corresponds to the prediction 
of the flat matter-dominated unperturbed Universe and the 
short-dashed purple line to a flat unperturbed Universe with 
fl A = 0.7 and Cl M = 0.3. 



the bare Hubble constant IKo with some physical prior. 
However, predicting a large variance for the deceleration 
parameter can be already regarded as an indication that 
the tree-level value go = \ ma y n °t be in conflict with ob- 
serving a locally accelerating Universe. A large variance 
does not imply the breakdown of perturbation theory ei- 
ther. Indeed, the density contrast Sp/p ~ ((pW 2 ip/a 2l X 2 ) 
(valid in the synchronous gauge for all scales) is at most 
of order unity. We have checked that a nonperturbative 
approach - along the lines of Ref. |l2( - leads to simi- 
lar conclusions and that yV 2 <£_is in fact the first term 
of the expansion of e~' p V 2 (p [l5|]. To get a physical in- 
tuition of the reason why there is a large contribution 
to the variance of the deceleration parameter one might 
think in terms of the energy density contrast. Indeed, 
a large variance of the <^V 2 y term implies a large vari- 
ance of Sp/p when summing up every Fourier mode of 
the perturbations. Such a variance is dominated by the 
infrared part of the spectrum and consequently is seen 
by an observer restricted to a region of size comparable 
to the present-day Hubble radius as a classical energy 
density background. Such interpretation is however not 
entirely correct because we are dealing with variances and 
not with averages of the physical observables. When in- 
specting the left-hand side of Einstein equations, one may 
regard the contribution from the first- and second-order 
gravitational potentials as kinetic energy of the gravita- 
tional field. Bringing them to the right-hand side and 
upon averaging, one may think of such terms as con- 



tributing to the classical energy density. However, this 
reasoning does not directly apply to variances. 

Instead of performing an angular average over the solid 
angle, we might have determined the expected values and 
variances of the Hubble rate and the deceleration param- 
eter as a function of the direction of observation. How- 
ever, since the variance is dominated by the infrared part 
of the spectrum, we expect that the local anisotropies of 
the variance are as small as {V 2 ip / a 2, K 2 ) . 

Finally, one should extend our computation to values 
of the physical redshift larger than unity. Indeed, an- 
other, albeit indirect, evidence of the presence of a Dark 
Energy fluid, is the measurement of the location of the 
first Doppler peak of the CMB anisotropies. Since the 
theoretical prediction for such a location depends, in an 
unperturbed Universe, on the total energy density pa- 
rameter J7o , from its measurement one may infer that Qq 
is very close to unity. Since Dark Matter amounts to only 
30% of the critical energy density, it is concluded that a 
consistent fraction of the energy density of the Universe 
is made of Dark Energy. However, this deduction holds 
in an unperturbed Universe only. On the contrary, one 
should compute the location of the first Doppler peak in 
a real perturbed Universe. This amounts to computing 
the angular diameter distance oIa subtended by the sound 
horizon on the last scattering surface. However, since the 
angular diameter distance is related to the luminosity dis- 
tance by the reciprocity relation cLa = + z) 2 at any 
order in perturbation theory, it suffices to determine the 
luminosity distance. The expression l|24|) for the luminos- 
ity distance-redshift relation holds for any redshift and 
any direction of observation and is therefore suitable to 
perform such a task. Furthermore, extending our compu- 
tation (for instance to higher redshift) is required to com- 
pare the theoretical predictions of a real perturbed Uni- 
verse to other physical observables supporting the Dark 
Energy picture in the unperturbed cosmology, such as 
the Integrated Sachs- Wolfe (ISW) effect and the transi- 
tion from the decelerating to the accelerating phase at 
redshifts of order unity 15]. In this respect we notice 
that the contribution to the variance of physical observ- 
ables, as V Var Kg)q l i ncreaS es with time and therefore 
both the ISW effect and the transition from the deceler- 
ating to the accelerating phase might well be consistent 
with a perturbed flat matter-dominated Universe. We 
also point out that the generic expression i|24|) may be 
used to estimate the influence of lensing on the brightness 
of supernovae sources. Finally, let us reiterate once more 
that a deeper understanding of the physical significance 
of long wavelength perturbations is certainly needed. 
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Appendix A: THE LUMINOSITY DISTANCE - REDSHIFT RELATION 

1. The geodesies equation 

To perform the calculation, it is useful to renormalize the photon four-momentum is such a way that 

g^k v (A = 0) = 1 . (Al) 

This is achieved by shifting the original photon momentum by a factor —1/ [ui(0)a(ri(0))], where r](0) is the value of the 
background conformal time at the observer point. . From now on, we will work with the shifted photon momentum. 

Expanding fc M = fc^s + fcj^ + k^ 2 y Eq. ijAlfl and the null normalization condition for k^ (the second of Eq. ©) 
imply the following initial conditions: 

fcf 0) = (-l,e% <%eV = l, k° {1) (0)=k ( { 2) (0) = 0. (A2) 
Expressing then ^, r — 1, 2 in terms of the perturbations of the metric, equation 1|22[1 gives 
dk?„ i . I 1 rl 2 



^(i) ~ \ (X(i)y)' eV = ^,y' eV (% - A) = \ (?/o - A) ^ , (A3) 



; + 2e* (e^ (1) ) - 2e l (1) )' - SkXfoeV + e fe ( X fo fc ) + -0* X(1)fcj . e V = --<p>* + y ^-eV 

I (Vo-^f ^ J C U 2f„ ^ 5„ i4 , 10^ , ( m -\) 2 d*(^) ^ 

dk?, 



dX 

fc 6 ^<, fc eV - - A)e^ = + fe<% + ^V^^ ~ - A)^ , (A4) 

°(2) 1 „,./ 1 



V»( 2 ) - J (x(2^)'eV + 2?A( 1) e l fc (1)i - (x(i)ij)' e t k {1)j + ^(i^fi) - i<9 M (^feV^j = 
(770 - A) p'V.i + (770 - A) 3 (<p M (p,ij - (VV) 2 ) - ^ ~ A) 3 ip'^^eV + ^ (770 - A) (p^e*)' 



rfA 2 
_5_ 

~ 18 

2 12 1 1 

+ (»7Q ~ A) V 2 ip^^e l e> - -S^eV + - (?/ - A) (p^e'k^ + -e fc eV,fcj?7(i) + ^eV^/yi^ (170 - A) = C(A) . 

(A5) 

In addition to the initial conditions for k9 r <., r = 1,2, given by equation l|A2|) . we can impose, first of all, 

afo (0) = xl 1} (0) = x\ 2) (0) = ?7(!) (0) = ry (2) (0) = , ry (0) (0) = m ; (A6) 

the initial conditions for fcl-g and fc? 2 ) , instead, must be chosen taking into account the null normalization condition. 
Imposing k^k a = for A = is indeed a necessary and sufficient condition to get a null geodesic, once k^.^k 1 ' = is 
solved [l3| . Assuming 

^ 1) (0) = ye l (A7) 
and imposing fc M (0)fc M (0) = to first order, we obtain 

*((>) (°) fc (o) (°) (- 2 ^(D ( x (o) M % + X(i)« (*(0) ,Vo))+ 2^o) (0)*J 1} (0)% - 2fc° 0) (0)k° {1) (0) = , (A8) 
from which we have 

1 -5 ri 2 ■ 

y = 0(i) (x( ),Vo) - ^X(i)»j («(o),»7o) eV = (x( )(0)) + ^e'e 3 ^ (aj(o)(Q)) ■ (A9) 
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We are now in a position to find the zeroth, the first and the second-order quantities necessary to define the photon 
trajectory 



V(o) = Va - A , 2( 0) = Ae 1 



= \ (<P(*(f>)W) - v(*(o}(0))) + km - A')eV,i(x (0 )(A')) 



A'=A 



A'=0 



(A10) 
(All) 



k\ t) (A) = ye* -2 /"^A' + ^eV(x (0 )(A')) + ± (fao - A') 2 ^ (* (0) (A')) e 3 - 2 V < i (x {0) (X'))(r }Q - A')) 

JO 6 A'=0 D V 7 



A'=A 



A A'=A 

»7(i) M = f j[ V ( x (o)( A 0) ^A' - ^ (x (0) (0)) - 3%AeV,i (sc«»(Q)) + ^ ((% - A') <p {x (0) {\'))) 



A'=A 
A'=0 

(A12) 
(A13) 



A'=0 



X (J X P* (*(o)( A ")) d\"\ dX' + |</(«(0))ijbA + ye* ( / r - (*„„ ( A' ) ) ,/A' - A r - (*„„ (())) 



z l (1) (A)=yAe l -2 
+ ~( - vlW 3 (x (0) (0)) ej + ((7? - A') 2 y>'' (x (0) (A'))) 



(A14) 



*° 2 )W = / C{X')dX' , 
Jo 

r K2)W=£ [J" C(X")dX") dX' 



(A15) 
(A16) 



Here and afterwards, if the argument of a function is left unspecified, we mean it is evaluated at (?J(o)(A), a;? )(A)), 

,2 

given by (|A10|) : for instance tp^je l e° = tp ^j(x^(X))e l e 3 = j${x(a)(X)). 



2. The optical scalars 

Let us now consider the transport equations for the optical scalars. To solve them, we obviously need to impose 
some boundary conditions. One could wonder if they are independent of the initial conditions we have imposed on 
the geodesic crossing the observer, Eqs. (jA2|>. IjAGf) . (|A7I) and (|A9I) (see Figure 0). Therefore, in order to understand 
which conditions must be adopted, let us suppose that our source emits isotropically (more precisely, suppose that 
an observer comoving with the source sees it emitting isotropically). Consider, then, a set (rj,x l ) of locally inertial 
coordinates around the source, associated with the tetrad 



{e (0 ), e ( i), e (2) , e (3) } , e (0) = u^-^ , 



(A17) 



where em), em, ep), e/3) are tangent to the space-time in the event corresponding to the center of the source; in 
particular, u^g|r = ^ is the 4- velocity of the center of the source. In other words, using tilded indices for the 

components in the coordinates (rj,x l ) and indicating with x c the coordinates of the center of the source, we require 
that, in a neighborhood of x c , 



9Jlv(. x ) = Vv-u + °(( x ~ x c) 2 ) and 



d 
dx^ 



Hp) 



(A18) 



(this yields, in particular, T~~(x) = 0(x — x c )). If we consider the geodesies on which we imposed IJA2I) . IJA6I) . (|A7|) 
and l|A9|l . in the center of the source the wave 4- vector reads 



with 



k^(X s + AX s )=u;(l,m l ) , 
u = -u^(X s + AA S ) = -1 + k° (1) {X s + AA S ) + fc^ 2) (A s + AX S ) : 



(A19) 
(A20) 
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Observer 




Figure 2: A congruence of photons emitted by a spherically symmetric source: one ray reaches the observer. 
m % eu\ (with m l m J dij = 1) is of course the direction in which the photon is emitted. 

If the emission is isotropic, the wave 4-vector of a photon emitted in a generic direction r^eu\ (with n l n^Sij = 1) is 
given by 

fc^(A s + AA s ,n) =^(1,^) . (A21) 

As we now have at our disposal the initial conditions (i.e. at A = A s + AA S ) for all the geodesies of the congruence, 
let us reconstruct it in a neighborhood of the source. First of all, let us note that 0(A — (A s + AA S )) = 0(2? — x c ): by 
definition, indeed, we have — kt 1 , which yields 

2" - 5£ = k»(\ s + AA.„n)j n _ »-» (A - (A s + AA S )) + 0((A - (A s + AA S )) 2 ) . (A22) 

Using therefore 

rg s (A) = 0(A-(A s + AA s )), (A23) 

it is immediate to get 

fc^A, n) = k»(\ s + AX s ,n) + 0(A - (A, + AA S )) 2 , (A24) 

which implies 

W - = k^(X s + AA„,n)| 5-a (A - (A, + AA S )) + 0((A - (A s + AA S )) 3 ) ; (A25) 

thus we have 

k*(x)=w(l,rj)\ n _ i-i +0((x-x r ) 2 ) . (A26) 

From this equation we easily obtain 

0{x) = Vjik* = 2UJ „ + 0(5 - xc) , (A27) 

\\X - X c \\ 

a(x) = 0(x-x c ) , (A28) 
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which, along a geodesic, become 

^(A)^ x - T ^^ y + 0(A-(A S+ AA S )), (A29) 

cr(A) = 0(A - (A s + AA a )) (A30) 

(use l|A21|) and l[A25)l ). It is therefore clear that the boundary conditions to impose on the transport equations for 
the optical scalars are 

^(o)(A)~ A _ 2 + as A^A S + AA S , (A31) 

a (0) (A s + AA s ) =0, (A32) 
a (1) {X s + AX S ) = 9 {1) (X S + AX S ) = 6 (2) {X S + AX S ) = . (A33) 

In particular, from (|A32|I we immediately get (see equation 114(1 ) 

tr (o) (A) = 0, VA; (A34) 

thus equation (|13fl becomes, to zeroth order, 



dX 

which yields, once Eq. I(A31|I is imposed, 



(A35) 



? (0) (A, AA S , A s ) = A A _ AA . (A36) 



Similarly, perturbing (|13|) to first and second order, we obtain 

^ = -0(o)0(i) - R(i)m»Wo) . (A37) 

from which we have, using Eq. 11331 . 

(1) (A, AA S , A,,) = — - _^ ^ (A' - A, - AA S ) 2 (%) Mi ,fcf 0) Afo)) (X')dX' , (A38) 



and 



~dX~ 



? ( o)0(2) (A, AA S ,A S ) 

' V (0) re (0) "I" ^- n -(l)^«( ) /t (l) f rl (l)A'^a x '(l) ft; (0) ft; (0) 1" 



S (A, AA S , A.) = R ( 2)^K k- a) + 2R (1)ll X ' fcfo + ^iSj^C + + 2|a (1) | 2 , (A39) 



with 



#(i)oo = 3^ ( i)" , (A40) 
>(i) = ^(x(i)t)' + 2 K)J'' ( A41 ) 

RiDij = %vV ( i) - ^v 2 X (i) y . + gxwjy + 2*MiM ~ 6i ^w" + + ' (A42) 



R 



(A - A s ) A 



2 



9(a) (A, A.) = - \ . 2 / (A'-A s ) z B(A',A s )dA' , (A43) 



B(X,X S ) = R(2)nvtf 0) kf 0) + 2i? (1)Aliy fc[ l 0) /c| y 1) + R^^xf^k^kfo + — + 2|cr (1) | 2 , 

11 3 

^(2)00 - 2 (x ( i) H )' (x ( i) H )' + M) kl (X(1) W )" + 3 (V(i)') + 6V(i)^(i)" + ^(2)" , 



I?, 



0*(2) 



m k i (^(D, fc ) - ^(i), fc (x(d ■ )' + V'(i) (x(i)l k ) - \x(D k 3 (xd)l fc ) 



+ \{*Wkj) X(i) fcJ ,i~^(i)* J 1 j(X(i)ifc)' + ^(i) fcj -(x(i) fcj ii ) +|( 
+ X(i)J fc ^(i)' + 4^(1) + 4^(i) (V(i) J + (V(2) J , 



X(2)i 



J2( 2 )« - % (W-(i)) 2 +X(i)}^(i) >w +X(i)-^(i) ifej +2%V(i)V 2 V'(i) 
- % + \ %VV(2) - |^(i) ,fc X(i)y ifc " ^(i)V 2 X(i )y . 

+ gX(i) X(i) tfifcI + 2^(i),fcX(i) i;j + ^(i)X(i) i;fcj ~ 2 X(1) ^ X ( 1 )ifc,b' 

_ ^X(i)J,X(i)5 ifc - ^X(i) fe 'X(i) j(ifei + \x(i) kl ti X(i) kU - *«^(i) ifc X(i) W iI 
+ 2*(i) (X(i) ijife - X(i) 4fej - X(i) jfcii J + 2 X(1) X^kl.tt ~ 4 V X ( 2 )*i 



1 



1 



+ \x(2) k M + \X{2)) M + 8 lJ (^(i)') 2 + X(i)l k i>(i),i + X(i) Jfc^i), 



+ 



1 1 1 / \ ' 1 / \" 

3 ^(i),i^(i) ,j + 2 ^(i)^(i), l3 " 2^^ (2) " + 2^ (2) .« + 2^ (1) ' V X(1) *J + 4 l X(2) *J 



^(i)oo(A) - -VV 



i? (2) oo(A) = -^V 2 (^ 2 ) + + ^V 2 * + ( Vo A) 2 (-^^°V,a6 + ^ (V 2 ^) 2 ) , 



i2(i)Oi(A)=i2(2)oi(A)=0, 
5 4 



55 50 /13 1 2 

^(2)«(A) = % y^'V.a + g-v 2 w*« + too - A) 2 IggP'V + yi ( V V) 

1 / \\4 / ,afcc ,ab , c , ,abc r w2 ,6c \ ^5 2 / 2\ 

^WO-A) (f,abc-<P' fi <P ■ ,ca + f,abf , c - %V W ' hc J - y V (<p) 

1 /a2 20 40 20 oo - 2N . 4 . 2 / 38 a 11 2 

- (a„ - V ) mj + —<p,np,j + -^-<P<P,ij - ~^d l d 1 [if ) + (no - A) ( -— if ti tp, aj - — V WA 



+ ^o-A) 4 



3 r^rj g rr, tJ g 

.,6 i .,,ab,„ i Y72,„,„,f> o .Ab 



63 



+ 21 -ip,aj<P' , bi -V,ai<P ,bj+ l P ,ab ( P^)+^ , b f \a: 
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Equation i|A39|) , using Eq. I|A33|) , is therefore solved by 



1 



A„+AA S 

2 



9 {2) (A, AA S , A s ) = -2 / (A' - A, - AA S ) B (A', AA S , A s ) dX' . (A53) 

(A-A S -AA S ) J a 

As B (A', AA S , A s ) depends on an) too, we must solve the first of Eq. 1)14(1. perturbed to first order, 

— ^ = -^(O)O-(I) + ^(X)ffM"« ft {0)*(0)*(0)*(Q) ' ( A54 ) 

1 / \" 

^(l)iojo = %^('i) - 2 V X(1 W ' ( ' A55 ^ 

(i)oyfe = fa (V'(i)j) - <*« (V»(i) jk ) + 5 - \ (x(i) iktj ) , (A56) 

1 / \" 



i? 



with 



ft 



ft 



2(i)o«fc = S ik - % (^(i), fc ) + 2 (*(i)<i,k) _ 2 ( X(1) ^j) ' (A58) 

«(i)<iki = &jrt{i), ik ~ Wfi),* ~ S ^W jk + S ^(i) s ~ \x{i) ik>jl + \x(i) lUk + \x{i) jk<a ~ \x{i) jhik , (A59) 
1 5 

R(l)iOjO = 3^>*J ' ^(l)Oij'A = ,R(l)ijkl = g (V,ifc^7 + fjtfik ~ <P,jk$il ~ f,il$jk) i (A60) 

with SL obeying, to zeroth order, the second of Eq. I|14fl : 



or 



"(0) 

u^ 1 + iw M 



= (0, «*) , tfyvV = 1, to" = (0, w l ) , JytoV = 1, 



^ V2 

r?7J l fill) 1 

S ljV l e 3 = S ijW z e J = <%iuV =0, — = — = . (A61) 

dX dX 

Imposing equation i|A33(l . the solution reads 

i /-A 3 +AA S 

ct (1) (A, AA S , A s ) = - — -—5 / (A' - A, - AA S ) 2 (^d^^o)^)^)^) ( A ') dX ' ■ (A62) 

(A — X s — /\X S ) J\ ^ ' 

3. The luminosity distance 

To get the relation between the physical radius of the source R and AA S let us use again a set (rj, x l ) of locally 

inertial coordinates in a neighborhood of the center of the source and let us define fc M = a~ 2 fc M (this corresponds to 
Eq. (J5J in the old affine parametrization). We therefore have 

M*) <*■> - (> - *&> < A '» - •&> < A '») " (I <*•>) (A63) 

(in particular, we have used Eq. JIJ and the new definition of /c M ). For the second of Eq. 10 yields drf 2 = Sijdx l dx^ , 
our relation reads 

R = ^5 l3 dxHxi = |A?7| = a {r, (A s )) AA S (l - fc° x) (A s ) - fc[> 2) (A s )) . (A64) 
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We can now solve for the amplitude A 



and find 



c?l = exp 



(Aa) (A.) = (Aa) (0) e'^o" e W dx ; 



(0 (o) (A, AA S , A s ) + 9 (1) (A, AA S , A s ) + (2) (A, AA S , A s )) dA 



(A65) 



x (1 + z(X s )) a AA s (l - fc° j (A s ) - fe° 2) (A,)) , 
where we have used a(ry(0)) = a(?jo) = ao- Taking the limit for AA S — > 0, noting that 



e -*/ o *»(0(o)(A,AA.,A,))dA = AA S 



and expanding the exponential, we get 
d L = (1 + z(X s )) a X 



AA S 



9 ( i)(A,A s )dA 



(A66) 



(A67) 



1 



9 (2) (A, Aa) dA - k° {1) (X s ) - fcp 3) (A.) + -k° {1) (X s ) / (A, A s ) dX 



, (A68) 



where 0m (A, A s ) and 0( 2 ) (A, A s ) stand for 0^) (A, AA S = 0, A s ) and 6r 2 ) (A, AA S = 0, A s ). In order to express the 
right-hand side of equation (|A68|) as a function of the redshift we can invert, order by order, the following equation 
(to get it, we have used Eq. JT5J|, (jA"2|) and JXgJ): 



l + z(A s ) 



Qq {u^) (X s ) ^ 
a( V (X s )) K^)(0) ~ 

ao(l (A s ) + fc° 2) (A.))) 



a (?7(o) (A s )) + a' (ry (0) (A s )) (?7 (1) (A s ) + rj {2) (A s )) + ±a" (?7 (0) (A s )) CV 
a 



a (V(0) (As)) 



(i + ri(A a ) + r 2 (A.)) 



(A69) 



(A70) 



t m I a '^(Q)) to a "0?(Q)) 2 , / a '(?7(o)) \ , a '(?7(0)) , \ m 



Adopting thus the definition 

A s = A( ) + A(i) + A( 2 ) 
and inserting it in the right-hand side of equation \ A69|) . we obtain 

a Q (1 + Ti (A (0) ) + T 2 (A (0) ) + T[ (A (0) ) A (1) ) 



l + z = 



A 2 , 



(»7o - A( )) - a' (770 - A (0) ) (A (1) + A (2) ) + a" - A (0 )) 



a (i70 - A (0) ) 



1 + a ' (^o-A(o))(A(i)+A (2) ) _ a" (yo - A (0) ) A^ ^ ^ 

a(»7o-A(o)) 2a(?7o-A(o)) 1 0) 



+ r 2 (A ( o))+T 1 / (A (0) )A (1) + 



a'(»7o-A(o)) \ a'(?7o-A(o)) 



a (% ~ A( )) 



Hi) 



a (r]o - A( )) 



A(i)Ti (A (0) ) 



(A71) 



(A72) 



(A73) 
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Setting then the zeroth order term equal to 1 + z and those of higher order to [5j , we have 



1 + 5 = 



(?7o-A ( o)) 



x (X ^_ Q (^0- A (0)) 

\i)V\a)) 77 : — 7 



T i ( A (o)) = r ?(i) ( A (0)J + 



- A( )) 



0.(77) 



(A74) 



(A75) 



A (2)(A( )) 



^ - A( )) 



T 2 (A (0) ) + T{ (A (0) ) A (1) + Ti (A (0) ) A (1) 5£ (% - A 



(o)J 



x U ^ A W a " (^o- A (Q)) 

+ [pt l^o - A (0 j A a) j --^-—7- 

^ a - A( ) J 



(A76) 



This procedure amounts to identifying the physical redshift parameter with the redshift in the unperturbed Universe 
evaluated at the zeroth affine distance Q 



z{X s ) = z(A (0) ) . 

If we replace A s in equation l)A68jl by means of l|A72(l and we expand in series, we obtain 



(A77) 



d L = 



a (r}o - A( )) 



(A( ) + A ( i) +A (2) ) 



A (°> 08 



1 r x (P) i 
1 ~ 2 J Q 9{1) ^ A(0) ' ) dX ~ 2 A(1) Jo UX 



(i) 



(A, A s ) |a s =a (0) ^A 



+ 8 U ^)( A ' A (°)) dA 



1 f X (0) 



) (2) (A, A (0) ) rfA - fc^j (A (0) ) - fc° 2) (A (0) ) - (A (0) ) A 



Mi) 



+ 2 fe (i)( A (0)) / ^(i)( A > A (o)H A 



d L (o) \ 1 



A(i) i r<°> 

A (o) 2 Jo 



hi) (A,A (0) )dA- fc^j (A ( o)) 



A( 2 ) 1 /' A(0) ##(1) ^(1) A(!) 1 / / ,A (°' ^ 

" ~ A (!) / Tn - ( A ' As ) U s =a ( o)^ a - -TT- (A( )) A(i) - - — k {1) (A (0) ) + - / 6» (1) (A, A (0) ) dA 



(o) 

1 / >A (0) 



(i 



9A S 



i) 



1 A f 



. (2) (A,A (0 ))rfA-fe[ , 2) (A ( o ) ) + ^ 1) (A ( o)) / " d {1) {X,X {Q) )dX-^ I " 6 {l) (A, A (0) ) dX 



(o) Jo 



The zeroth order luminosity distance is given by 



d H0) = 7 a ° s A ( o) = — (l + z - Vl + z) 
a (//o - A(o)j "H-o v / 



(A78) 



(A79) 



where we have used a (77) = a ( ^ ) and the inverse of Eq. (|A74|) , 



A (Q) _ 1 1 
Vo 



y/T+z 



(A80) 



assumption that R^^k^k^, B (X' , AA S ,A S ) and R^a^sk^k^t^t^ can be written as a power-series 



Besides, it's worth noting that, in order to get equation l|A78|) . we have used 6'( 1 )(A( ), A( )) = 0. Making indeed the 

+00 

^c„(A-A s ) n , (A81) 



n=0 
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2(1); 0(2) an d take the form 

n + 3 



n=0 



Appendix B: THE DECELERATION PARAMETER 

As already anticipated, to get % we can replace \o)/vo = 1 ~ l/\/T+~z in Eq. (|A78|) and expand in series 
around z = (for simplicity, we will set % = 1)- Denning 0(i) (O,1) (O, 0) = |^(A, A s )| a =a s =o and 6> (1) (1 ' 0) (0, 0) = 
gf(A, A s )| a=a s =o an d indicating with a prime the differentiation of rjm and 7j( 2 ) with respect to A (for example 
J7( X \(0) = -^HA)|a=o = /c^(0)) and that of A(i) and A( 2 ) with respect to A(o) (see equations (|A75|) and (|A76|> ). we 
obtain 



A + T,I + 6? (Bl) 



with 



A = A (1) (0) + A (2) (0) - A (1) (0)^(0) , (B2) 

3 = | + A (1) (0) - + A(2)(0) _ _ A(1)(0) ^ )(0) 

^( 2 )(0) A' (1) (0) ^ 1} (0)A' (1) (0) A' (2) (0) 



2 2 2 2 



A(i)(0) V(i)(0), (B3) 

e = i (i - 0a) (o, o) - e {2) (o,o)- 0(1)( °'° 2 )A(1)(O) - (o) + e {l) (o, o) v ' (1) (o) 
-»?(3) (°) + A (D (°) " e (D (°> °) A (D (°) - ii) (°) A (D (°) + A w (°) ~ 2 <d (°) 

-2A (1) (0) t/fa (0) - 4A' (1) (0) t/fa (0) - 2^' 2) (0) + A' ( ' x) (0) - (0) X'fa (0) + A' ( ' 2) (0) 

-3^(0)^(0) - 2A (1) (O)0 (1) (O ' 1) (O,O) - ( H (1 j --J . (B4) 

Using equations l)A75(l and (|A76|) together with 

(0) = V(i) (0) = r? (2) (0) = »7(2) (0) = 0(i) (0, 0) = (2) (0, 0) = (B5) 

(the first four relations are among the initial conditions we have imposed on the geodesies equation, while the last 
two follow from equation (|A82(l 'l. we can easily get 

A = 0, (B6) 
3 = i 



I ^4 + 2rf(y ) (0) + (y ( ' 1} (0)) 2 + 2r?' ( ' 2) (0)) , (B7) 

I (4 - 6 r/fa (0) + 2 ^ (0) - 4 r/ ( ' 1} (0) 2 - 6 ^ (0) 
+3W'(0)<i')(0) + 2< 2 ' ) (0)) . (B8) 



e 

32 



Comparison between Eqs. (|51Jl and (|B1|I then yields 
1 

qo =8 



4 + 2 (4^ (0) - ^ (0)J + < x) (0)" + 8,/ ( ' 2) (0) - 2^ (0)^ (0) - 2r/fc (0) j , (B9) 
JC = 2-7 ? ( / 1) (0)-ry| 2) (0). (B10) 
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We can express all the parameters in terms of the original metric 7^ as 9 

<i)(0) = -^(7(i W )'eV, (Bll) 

<i')(0) = \ (7(i)«)" eV \ (0 fc7(1)ij .)' eVe* , (B12) 
< 2) (0) = -| (7(2)«)'eV- - ( 7 (i W )VfcW(0) , (B13) 
<2)(0) = 1 (7(2)«)"eV- - i (0 fc 7(2)«)'eVe fc + (7(i)«) VfcUo) 



^ fc (^DoO'e^uW - ( 7(1)y )V^i(0) - ^ ( 7(1)ij )' eV^^O) , (B14) 



^(0) = -d kl \ 1)3 e k e> + e k ( 7( \ )fe )' + ^ 7(1)jy eV , (B15) 

4 ) (0) = -^ 7(1)fej eVe i . (B16) 

The functions 7 ( r ) (r = 1,2) and their derivatives are evaluated at the observer's position X( ) = at the present 
time r]o = 1. In terms of the gravitational potentials we have 

<i)(0)=V-(i)-^(X(i)«)'eV, (B17) 

<i')(o) = -V('d + ^ (xdw )" eV ' + e*Mu - ^ fe eV ' - ( B18 ) 

<2)(0) = \^2) - \ (X(2 )i ,)'e i e j + 2^ 1) e i fc (1)i (0) - ( X(1)tf )' ^^(0) , (B19) 

<2)(°) = -^(2) + J (X(2 W )"eV + ie fe a fc ^ 2) - \e k (d fe x (2)jj )'eV - 2^' 1) e i fc (1)i (0) + ( X(1) «)" ^(O) 
+ 2e^ 1) e i fc (1)i (0) - e k d k { X[m )' ^(0) + 2^/^(0) - e'^(0) + ^^,^,(0) 

-^(^Xd^O'^e^W, (B20) 

,(1) -(0) = + 2e* (e^ (1) ) - 2e* ty (1) )' - SkX^-eV + e fc (V (1)fe ) + ^Xd^eV , (B21) 



fc(i)(0)= ( V-(i) - Jx ( i) fej eV ) e\ (B22) 



9 



When applied to the metric, the prime stands for differentiation 
with respect to conformal time: ' = d v . 
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where ipi r \, X(r) ( r = 1, 2) and their derivatives are evaluated in the observer's position imi = at the present time 
Vo = !■ 

What we need to make a comparison with the observations are the averages of these quantities over the sky (that 
is, over the direction of observation). 

Let us call m the direction of observation (by definition, m • m = 1): we can of course expand m either in the 
coordinate basis {gfrj g§r? g§s} or m a Cartesian system of axes ^(2)i n (3)} ( a "triad" carried by the observer): 

m = ™^i=^W- (B23) 
In our case, the e 1 , i = 1, 2, 3, are proportional to the m l , i = 1, 2, 3: 

e i =Nm i , (B24) 



where the constant N is determined by the normalization condition Eq. (|A2|) . In order to perform an angular average, 
we therefore need to express equations l|B17|) - l|B22l) not in terms of the e l , but in terms of the director cosines e\ 
Let us therefore define 

d 1 d 

n {k ) - (1 + V>(i))^fc ~ 2 X « fc fo? ; (B25) 

to first order we have 

n (i) ' n ti) = % : ( B26 ) 

the Tim, i = l,2, 3, are the unit vectors of a triad. 

Taking then into account Eqs. 1|A2|). [jB23j> . (|B24|) and (|B25|) we get, to first order, 



e l = V 'gkie k e l = \J gue k e l m l 



1 kl 



1 - -0(1) + ^X(i)efeei 



'''' ' " ^ (xf^fcQ) • (B27) 



Performing this substitution in equations (|B17|) - (|B22I) . we can calculate the angular averages of 3€o and qo by means 
of the following identities (which can be obtained by expressing the director cosines e* in spherical coordinates): 

<e% = (eVe*> n = (eV^)o - , (B28) 

2n+l terms 

(eV) n = , (B29) 



(e i e j e k e l ) u = — (6 lj S kl + 5 lk S jl + S ll 5 k: >) , 



(B30) 



where we have indicated with ( )q the average over solid angle 



-in 



(B31) 



Therefore, by using these identities and equations (|A3(I - (|A5(I . the angular averages of equations l|B9|l and (|B1 

become, after resurrecting the dependence on t}q, 



1 " Vo QV(i) + 1^(2) + ^(i)V'(i) + ^ (X(i)) X(i)i 

ivv - ^ (v,)' + ,S - + ^ (vV) J ) ,5 



(B32) 
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(B33) 



where tp( r ), X(r) i r = 1, 2) and <p are evaluated at the observer's position CE( ) = and at the present time 770 and 
where "Kq an d 9o are the background Hubble constant and deceleration parameter IHo = and q = \ . 
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